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Background: In the search for evidence of life, 
Icebreaker[l] will drill in to the Martian ice-rich rego- 
lith to collect samples, which will then be analyzed by 
an array of instruments designed to identify biomark- 
ers. In addition, drilling into the subsurface will pro- 
vide the opportunity to assess the vertical distribution 
of ice to a depth of 1 meter. The purpose of this par- 
ticular project was to understand the uncertainties in- 
volved in the use of the imaging system to constrain 
the water ice content in regolith samples. 

Mixture Modeling. The key to constraining the wa- 
ter ice content of a mixture is to model said mixture. 
Shkuratov et al. modeling [2] was used instead of 
Hapke Modeling [3] [4] due to the eschewing of view- 
ing geometry in calculating the albedo (which can be 
related to the bidirectional reflectance[5]), as well as 
the fact that the Shkuratov model is reversible in that 
the imaginary index of refraction of a material can be 
calculated from its albedo. 

Subtractive Kramers -Kronig analysis. Knowing the 
imaginary index of refraction is of little use without the 
real index of refraction as well. By using Subtractive 
Kramers-Kronig analysis [6], one can use the real index 
of refraction at a specific wavelength to iterate through 
and derive the real and imaginary indices of refraction 
for the rest of the wavelengths in the range ob- 
served^]. 

Coming together. The idea then becomes that one 
can determine the fraction of water ice in a sample by 
taking a spectrum of a regolith sample, wait for the 
sample to dessicate, and then taking another spectrum 
using the same lighting and viewing geometry. From 
the second spectrum, the optical constants for the ice- 
less mixture as a whole can be determined using the 
Subtractive Kramers-Kronig analysis. These optical 
constants can then be used in conjunction with the op- 
tical constants for water[8] to model mixtures of icy 
regolith until it matches the spectrum of the observed 
icy regolith. Theoretically, this would determine the 
water ice abundance of the regolith, which in turn 
gives insight to the chemical and geological history of 
the region. 

A Case Study: Early experiments wherein samples 
were observed before and after sublimation of water 
ice seemed to support the theory that the water ice con- 
tent of these samples could be estimated using Shkura- 
tov modeling and Subtractive Kramers-Kronig Analy- 
sis. An example of this is seen in Figure 1. 


Modeled Spectra assuming n = 1.4 at 550 nm, 10 micron ice grains 



Figure 1: Silicate Sand was observed with two concentrations of ice 
mixed in. Spectra were taken of these icy samples as well as the 
dessicated sample. Above are the modeled mixtures overlaid on the 
spectra of the original samples. Grain size of the silicate sand was 
assumed to be 100 pm, with a volume fraction filled by particles of 
0.85. It was also assumed the silicate sand had a real index of refrac- 
tion of 1 .4 at 550 nm. 

For another test of the theory, Phoenix data was 
looked at, of samples before and after sublimation on 
the Martian surface — in particular of Snow White 
Trench between Sol 45 and 50 of Phoenix’s mis- 
sion^]. 

Unlike with the silicate sand observed in Figure 1, 
the viewing geometry of the Snow White Trench was 
not quite of incidence angle 30°, emission angle 0°; 
however, the formula log(7? 30 ) = 1.0 8 8 * log(^4) was 
used to approximate the hemispheric albedo used by 
Shkuratov et al. modeling from the bidirectional reflec- 
tance [5]. This approximation should remain accurate 
within a few percent reflectance. 

Unlike previous models that attempted to deter- 
mine the water ice content of the trench by assuming 
the composition of the regolith[10], we derived the 
optical constants from the observations of the sublima- 
tion lag using the Subtractive Kramers-Kronig analy- 
sis outlined earlier. Since the true real index of refrac- 
tion of the regolith is not known, it is assumed to be 
between 1.4 and 1.8 at 550 nm. The optical constants 
of water ice come from [8]. Various concentrations of 
water ice and regolith are then modeled using [2]. 

The volume fraction filled by particles was as- 
sumed to be 0.85. Varying this fraction does not shift 
the modeled spectra to a great degree — perhaps by 0.1 
albedo when varied from 0 to 99% porosity. The effec- 
tive grain size of the sublimated regolith was assumed 
to be 60 pm[10]. However, it does not matter what 
grain size is used to derive the optical constants of the 


regolith so long as this same grain size is used in the 
modeling of this regolith. 

Results: Various mixture models were generated, 
assuming different real indices of refraction for the 
regolith, as well as different ice grain sizes. Figures 2-5 
show these models alongside the original spectra (both 
initial and sublimated). 


Modeled Spectra assuming n=1.4 at 550 nm, 10 micron ice grains 



Figure 2: Modeled mixtures overlaid on actual spectra of the Snow 
White Trench before and after sublimation. Assumed n=1.4 at 550 
nm, regolith grain size of 60 pm, porosity of 0.15, and ice grain size 


of 10 pm. Suggests water ice content lies between 5-20% ice. 
Modeled Spectra assuming n=1.4 at 550 nm, 100 micron ice grains 



Figure 3: Modeled mixtures overlaid on actual spectra of the Snow 
White Trench before and after sublimation. Assumed n=1.4 at 550 
nm, regolith grain size of 60 pm, porosity of 0.15, and ice grain size 


of 100 pm. Suggests water ice content lies between 5-20% ice. 
Modeled Spectra assuming n = 1.8 at 550 nm, 10 micron ice grains 



Figure 4: Modeled mixtures overlaid on actual spectra of the Snow 
White Trench before and after sublimation. Assumed n=1.8 at 550 
nm, regolith grain size of 60 pm, porosity of 0.15, and ice grain size 
of 10 pm. Suggests water ice content lies between 10-30% ice. 


Modeled Spectra assuming n=1.8 at 550 nm, 100 micron ice grains 



Figure 5: Modeled mixtures overlaid on actual spectra of the Snow 
White Trench before and after sublimation. Assumed n=1.8 at 550 
nm, regolith grain size of 60 pm, porosity of 0.15, and ice grain size 
of 100 pm. Suggests water ice content lies between 10-30% ice. 

The main effect of increasing the assumed real in- 
dex of refraction of the regolith is to decrease the al- 
bedo of the modeled mixtures. Increasing the grain size 
of the ice grain meanwhile appears to have little effect 
for reasonable size ranges (however far larger in- 
creases would begin to lower the albedo in the near 
infrared wavelengths). 

The plots show that this water-ice content is likely 
between 5 and 30%, but the models are not perfect. 
Not knowing the real index of refraction (as well as the 
approximation of the hemispheric albedo from the 
bidirectional reflectance) creates much uncertainty in 
the water-ice content of the sample. This can be 
mitigated by a closer approximation of the hemispheric 
albedo (or more constraining the viewing geometry), 
and knowledge of the real index of refraction (which 
can be approximated if one knows the constituents of 
the regolith). 
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